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Introduction 


Comparison  with  measurements  is  final  goal 

Validation  challenge:  many  parameters  and 
assumptions  needed  to  model  shock  experiments 

Need  to  examine  them  for  chemistry  and  vibrational 
relaxation  models 

Recent  years:  many  detailed  QCT  computations, 
especially  for  N2-N2  and  N2-N 

These  data  provide  basis  for  benchmarking  and 
accuracy  analysis  of  fast  empirical  models 
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Objective 

Our  project  seeks  to  recommend  and  improve 
simplified  microscopic  models  and  connect  them  to 
macroscopic  counterparts 

This  work  is  first  step:  use  recent  high  quality  QCT 
calculations  for  model  verification  and  parameter 
adjustment 

Four  chemistry  models:  total  collision  energy 
(TCE),  quantum  kinetic  (QK),  vibration-dissociation 
favoring  (VFD),  and  weak  vibrational  bias  (Bias) 
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Numerical  Approach 

DSMC  method:  SMILE  tool  extended  to  include 
Bias  and  QK  models 

Homogeneous  bath:  N2  and  N 

Translationally  isothermal:  velocities  re-sampled 
every  time  step 

Majorant  frequency  scheme:  modified  to  correct 
for  larger  than  1  reaction  probability  in  Bias  model 

VHS  model:  d  and  co  to  match  high-T  viscosity 

Diameter  of  N2  depends  on  internal  state  (v,J) 

Internal  energy  modes:  continuum  rotational, 
discrete  vibrational  (AHO) 
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Distribution  Function 


Vibrational  Energy  Transfer 


VT:  Larsen-Borgnakke  model  with  ZV(T)  prohibiting  double 
relaxation  -  mathcing  QCT  for  all  T 
VV:  QC  model 


VV  reduces  unphysical 
overrepresentation  of 
ground  level  (transient 
state  at  TV=4,500K) 
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Reasonable  agreement 
with  QCT  (QSS  state  at 
T=30,000K) 
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Non-Coupled  Dissociation  Models 

TCE  model:  reaction  occurs  when  E=E+E  +E  >D 

try 

with  probability  ( e-dY 1 

p=di  - - — 
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ap  bp  c1  depend  on  Arrhenius  and  VHSA/SS 
parameters  and  molecular  properties 
TCE  is  tied  to  VHSA/SS  interaction  model 


QK  model:  reaction  occurs  during  VT  energy 
transfer  when 
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QK  is  integrated  with  Larsen-Borgnakke  VT  transfer 
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Models  Coupled  with  Vibrations 

VFD  model:  reaction  occurs  when  E-E  >D  with 

v 


a2,  b2,  c2  depend  on  Arrhenius  and  VHS/VSS 

<|>  is  free  parameter  that  controls  vibration  coupling 


Bias  model:  reaction  occurs  when  E.+E  >D  with 
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X  controls  vibration  coupling 
A  is  adjusted  to  match  thermal  reaction  rate 
Simplest  to  implement,  not  tied  to  any  other  model 
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Dissociating-to-Total  Vibrational 
Distribution  Ratio,  N,-N  at  QSS 


General  trend:  reaction  rate  increases  with  v 

TCE,  QK:  lack  of  vibrational  favoring  results  in  much  lower 

slope  as  compared  to  the  benchmark  QCT 

VFD:  (|)=1  only  works  for  lower  levels,  <j)=3,  for  high  levels 
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Vibrational  Distribution  for  Bias  Model 

10,000  K  30,000  K 


•  Excellent  fit 

•  10,000  K:  X=4  and  X=5  work  well 

•  30,000  K:  X=4  provides  best  fit 

•  Impact  of  diameter  increasing  with  v  and  J  is  small 
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Reaction  Rate,  cm3/s 


Total  Equilibrium  Reaction  Rates 


Adapted  from  Jaffe  et  al  AIAA  2016-0503 


•  TCE  and  VFD:  explicitly  match  total  reaction  rate 
written  in  Arrhenius  form 

•  Bias  model:  A  needs  to  be  adjusted,  but  captures  rates 
well  for^4~1 

Distribution  A:  Approved  for  Public  Release,  Distribution  Unlimited  PA  #16299 


Reaction  Rate,  cm3/s 


Level  Specific  N2-N  Reaction  Rates 


•  TCE,  QK:  lack  of  vibrational  favoring  results  in  much  weaker 
dependence  of  the  reaction  rate  on  v  as  compared  to  QCT 

•  VFD:  cp=3  provides  reasonable  agreement  with  QCT  for  v>30, 
but  for  the  first  few  levels  it  is  orders  of  magnitude  lower 

•  Bias:  very  good  agreement  for  X=2 
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Dissociation  Constant,  m3/s 


Non-Equilibrium  Reaction  Rates 

N2-N2,  T/T=Tr 


•  TCE,  QK,  and  VFD  are  over  an  order  of  magnitude  off  at 
lower  Tv 

•  Of  these  three  models,  TCE  provides  somewhat  better 
agreement  with  QCT 
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Non-Equilibrium  Rates  for  Bias  Model 

N2-N2>  T/T=Tr 


•  Bias  model  provides  much  better  fit  than  Park  2-T  model 

•  Higher  vibrational  favoring  for  lower  Tt:  X>2  works  better 
than  A=4  for  T>20,000  K 
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Reaction  Rate,  m3/s 


Non-Equilibrium  N2-N2  Reaction 
Rates  for  Low  Internal  Temperatures 


1  1.5  2  0.5  1 

10000/T,  1/K  10000/T,  1/K 


•  TCE  within  one  order  of  magnitude 

•  VFD  model  fails  at  low  Tv 
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Reaction  Rate,  m3/s 


Non-Equilibrium  Rates  for  Low 
Internal  Temperatures:  Bias  Model 
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1=4  for  N2-N2  and  X=2  for  N2-N  provide  good  agreement  with  QCT 
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Relaxation  from  Equilibrium  to  QSS 

Initial  state:  Tv=Tr=Tt=30,000  K 


Bias  model  agrees  well  with  QSS 
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N,  Concentration 


Relaxation  from  Non-Equilibrium  to  QSS 

Initial  state:  TV=T=2,000  K,  T=30,000  K 


•  QSS  state  for  1=2  is  close  to  QSS,  but  initial  rotational 
relaxation  is  faster  and  vibrational,  slower 

•  Possible  reason:  LB  model,  where  relaxation  numbers  are 
determined  from  local 

•  VT  rate  reduced  by  2  gives  better  agreement  over  the  first  ns 
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N,  concentration 


Conclusions 


Accuracy  of  four  DSMC  dissociation  models  is 
analyzed  through  comparison  of  N2-N2  and  N2-N 
with  recent  QCT  calculations 

TCE,  QK,  and  VFD  models  do  not  capture  key 
reaction  features  related  to  the  vibration- 
dissociation  coupling 

Bias  model  performs  very  well,  and  the  adjustable 
parameters  X  =  4,  A  =  1 .5  for  N2-N2  and  X  =  2,  A  = 
1 .5  for  N2-N  provide  good  agreement  with 

These  preliminary  results  set  the  stage  for  more 
extensive  comparisons  of  this  promising  reaction 
model,  to  include  more  realistic  cases 
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